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Abstract 

Given an undirected graph with edge costs and node weights, the minimum bisection problem 
asks for a partition of the nodes into two parts of equal weight such that the sum of edge costs 
between the parts is minimized. We give a polynomial time bicriteria approximation scheme for 
bisection on planar graphs. 

Specifically, let W be the total weight of all nodes in a planar graph G. For any constant 
5 > 0, our algorithm outputs a bipartition of the nodes such that each part weighs at most 
W12-\-e and the total cost of edges crossing the partition is at most (1 + ^) times the total cost of 
the optimal bisection. The previously best known approximation for planar minimum bisection, 
even with unit node weights, was 0(log n). Our algorithm actually solves a more general problem 
where the input may include a target weight for the smaller side of the bipartition. 
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1 Introduction 

Breaking up is hard to do. The most famous hard graph-breaking problem is graph bisection: 
partitioning the vertices of a graph into two equal-size subsets so as to minimize the number of edges 
between the subsets. This problem was proved NP-hard by Garey, Johnson, and Stockmeyer [12] 
in 1976. 


Background But how hard is it really? In particular, how well can graph bisection be approxi¬ 
mated in polynomial time? Even assuming P ^ NP^ we cannot at this point rule out the existence 
of a polynomial-time approximation scheme (PTAS) for graph bisection. The best approximation 
ratio known to be achievable in polynomial time is O(logn) where n is the number of vertices, due 
to Racke [24], improving on a bound of O(log^'^n) due to Feige and Krauthgamer [10] (the first 
result discovered that had a polylogarithmic approximation ratio). 

One way to make the problem easier is to relax the balance condition. Given a number 0 < 
b < 1/2, a bipartition U UV of a graph’s vertices is b-balanced} if \U\ > [bn] and \V\ > [bn]. Any 
bipartition U UV induces a cuh namely the set of edges between U and V. The bisection problem 
is to find a minimum ^-balanced cut. It might be a simpler problem to find a nearly optimal 
6-balanced cut for some b < 1/2. No better approximation ratio is known for this problem when 
the input graph is arbitrary. However, for the special case of planar graphs, a 2-approximation 
algorithm was given by Garg, Saran, and Vazirani [13] for finding a minimum 6-balanced cut for 
any b < 1/3. 

Bicriteria approximation^ gives another way to relax the balance condition. A bicriteria approx¬ 
imation algorithm seeks a ^'-balanced cut whose size is at most some factor times the minimum size 
of a 6-balanced cut. In an early and very influential paper, Leighton and Rao [19] showed, using 
a reduction to their 0(logn)-approximation algorithm for another problem, uniform sparsest cuf 
that a ^'-balanced cut could be found whose size is 0(^^) times the minimum size of a 6-balanced 
cut for any b' <b and b' < 1/3. Using the improved 0(\/log n)-approximation algorithm of Arora, 
Rao, and Vazirani [2] for uniform sparsest cut, the approximation ratio 0(^^) can be improved 

to Note that this performance ratio gets worse as the graph size grows and gets worse 

as the balance b' achieved by the algorithm approaches the balance b that defines the optimum. 


Our results In this paper, we give a bicriteria approximation scheme for bisection in planar 
graphs: 

Theorem 1.1. For any e > 0, there is a polynomial-time algorithm that, given a planar graph G, 
returns a ^ — e-balanced cut whose size is at most 1 + e times the optimum bisection size. 

That is, the algorithm returns a partition of the vertex set that is almost perfectly balanced 
(each side has at most a fraction ^ -h e of the vertices) and whose size is almost as small as the 
smallest bisection. 

Previously no approximation algorithm was known that had a 1 -\- e approximation ratio even 
if the algorithm was allowed to return a ^'-balanced cut for some constant F > 0, even if F was 
allowed to depend on e, even for planar graphs. The algorithm generalizes to handle 6-balanced 
cuts: 


^Some papers use this term to mean that each part has cardinality at most bn. 
^This approach is also called a pseudo-approximation. 
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Theorem 1.2. For any e > 0, there is a polynomial time algorithm that, given a planar graph G 
and given 0 < b < ^, returns a [h — e)-balanced cut whose size is at most 1 + 6 times the minimum 
b-balanced cut. 

Also, the algorithm can handle nonnegative costs on edges and nonnegative weights on vertices. 
In fact, we prove a more powerful theorem. Let G be a graph with vertex-weights and edge-costs. 
For a number b, a b-bipartition of G is a bipartition U UV of the vertices of G such that the total 
weight of U is exactly b times the total weight of G. The cost of the bipartition UUV is the cost of 
the corresponding cut, i.e. the set of edges connecting U and V. For example, a minimum bisection 
is a minimum-cost ^-bipartition where all costs and weights are 1. Let OPT 5 (G) be the cost of an 
optimal 6 -bipartition. 

Theorem 1.3. For any e > 0, there is a polynomial-time algorithm that, given b > 0 and given 
a planar graph G with edge costs and vertex weights such that G has a b-bipartition, returns a 
b'-bipartition whose cost is at most (1 + 6 )OPT 5 (G), where 6 ' G [b — e,b -\- e]. 

1.1 Related work 

There is much prior work on finding approximately optimal separators in planar graphs. Before 
the work of Leighton and Rao [19], Rao [25,26] gave approximation algorithms for balanced cuts 
in planar graphs. One is a true approximation algorithm for finding a 6 -balanced cut (for b < ^) 
whose performance guarantee is logarithmic. Another is a bicriteria approximation algorithm whose 
performance guarantee grows as the balance b' achieved by the algorithm approaches the balance 
b that defines the optimum. In our algorithm, we employ a subroutine of Rao [26]. 

Park and Phillips [23] gave improved algorithms for achieving some of the goals of [25,26]. 
The aforementioned result of Saran, Garg, and Vazirani [13], the 2-approximation algorithm for 
6 -balanced cut in planar graphs for 6 < 1/3, built on the work of Park and Phillips. 

There has been work on approximation schemes for other graph classes. Arora, Karger, and 
Karpinski [1] gave an approximation scheme for bisection in unit-edge-cost dense graphs, graphs 
with edges. Guruswami, Makarychev, Raghavendra, Steurer and Zhou [14] gave an algorithm 

that, given a graph in which there is a bisection that cuts a fraction 1 — 6 of the edges, finds a 
bisection that cuts a fraction 1 — g{e) of the edges, where g{e) = 0{^\og{l/e)). 

There has been much work on using approximation algorithms for balanced separators to obtain 
approximation algorithms for other problems; see the survey of Shmoys [28]. 

There has been much work on finding balanced separators of size 0{^/n) in planar graphs, 
regardless of the optimum bisection size. Lipton and Tarjan [20] first gave such an algorithm. 
Many papers built on this result: improvements to the multiplicative constant, an algorithm to find 
separators that are also simple cycles [ 21 ], algorithms that build on geometric embeddings [ 22 ] or 
eigenvectors [29]. In addition, there has been much work on algorithms that use planar separators. 


2 Overview 

We outline the algorithm for Theorem 1.3, the bicriteria approximation scheme for min 6-bipartition. 
The input is a graph G* with vertex weights and edge costs. The algorithm involves edge contrac¬ 
tion. Contracting an edge uv means removing the edge and replacing its endpoints u and v with a 
single vertex x. Edges previously incident to or 'L’ are now incident to x. The weight assigned to 
X is the sum of the weights of u and v. 
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2.1 Framework for approximation scheme 

The algorithm uses a framework of Klein [18]. The framework has previously been used to address 
optimization problems in planar graphs [4-8,18], such as traveling salesman and Steiner tree, that 
involve minimizing the cost of a set of edges subject to connectivity constraints. The framework 
has never been used before in the context of weight constraints on the vertices. 

First we give the outline. Fix 0 < 6 < 1 and e > 0. The framework uses the notion of 
branchwidth^ [18,27]. 

1. Spanner step: In the input graph G*, contract a selected set of edges, obtaining a graph G* 
with the following properties: 

• c(^) < pOPT 5 (G*), and 

• there exists 6' G [b — e,b -\- e] such that OPT 5 /(G*) < (1 + e)OPT 5 (G*). 
where p is a quantity that depends on the construction. 


2. Thinning step: Select a set S of edges in G* such that: 

• c{S) < (l/A:)c(G*) and 

• G* — 5 has branchwidth 0{k), 

where k = e~^p. 

3. Dynamic-programming step: Find a cheapest 6'-bipartition {Ui,U 2 ) in G* — S, where b' G 

[6 — e, 6 + e]. 

4. Lifting step: Return [Ui, U 2 ) where Ui is the set of vertices of G* coalesced to form vertices 
in Ui. 


The cost of the returned solution is at most 


OPT5/(G*-5) + c(5) < 

< 

< 

< 


OPT5/(G*) + c(5) 

(1 + 6)OPT5(G*) + (1/A:)c(^) 
il + e)OPTb{G*) + ep-^ciG'*) 
(1 + e)OPTbiG*) + eOPTbiG*) 


which shows that the cost is at most 1 + 2e times the cost of an optimal 6-bipartition. The fact 
that {Ui, U 2 ) is a 6'-bipartition of G* means that {Ui, U 2 ) is a 6'-bipartition of G*. 

The thinning step is straightforward: choose a root, and find breadth-first-search levels for all 
edges in G*. For i = 0,1, 2,..., A: — 1, let Si be the set of edges whose levels are congruent mod k 
to i. For each i, G* — Si has branchwidth 0{k) (see, e.g., [18], also treewidth 0{k), see, e.g., [3]) 
and at least one of the sets Si has cost at most (l/A:)c(G*). 

The fact that G* — S has branchwidth 0{k) means that in the dynamic-programming step an 
optimal 6'-bipartition can be found in time 2^(^)poly(n, W) where W is the sum of weights. 

Assume for now that W is 0{n) and that p is O(logn). Then the dynamic-programming step 
requires only polynomial time. (This is a straightforward generalization of Theorem 4.2 of [16].) 

The one challenging step is the spanner step.^ The main work of this paper is showing that this 
step can be done. 

^ Treewidth could be used instead 

^This is usually the case in applications of the framework. 
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Theorem 2.1. There is a constant c and a polynomial-time algorithm that, given e > 0, b > 0 
and a planar embedded graph G* with vertex weights and edge-costs, returns a graph G*, obtained 
from G* by contracting edges, with the following properties: c(G*) < clogn • OPT 5 (G*), and 
3b' e[b-e,b + e] such that OPT 5 /(^) < (1 + 6 )OPT 5 (G*). 

Once we have proved Theorem 2 . 1 , showing that there is a poly-time algorithm for the spanner 
step of the framework, we will have proved Theorem 1.3. 

2.2 Spanner construction overview 

Many tools have been developed for spanner construction. One tool first used for Steiner TSP is 
this boundary-to-boundary spanner: 

Lemma 2.2 (Theorem 6.1 of [17]). Let G be a planar embedded graph with edge-costs, and 
let G be the boundary of some face of G. For any 6 > 0, there is a subgraph H of cost 0(e“^c(G)) 
such that, for any two vertices x and y of G, the x-to-y distance in H is at most 1 + 6 times the 
x-to-y distance in G. Furthermore, there is an O(nlogn) time algorithm to derive H from G. 

The edges defining the min-cost bisection or 6 -bipartition of input graph G* correspond in the 
planar dual G to a collection of edge-disjoint cycles. Fragments of these cycles are paths; perturbing 
the solution by replacing such a path with a nearly shortest path in a subgraph H does not increase 
the cost of the solution by much. This simple idea is at the heart of the spanner construction. 

At the highest level, we use the following strategy: (Step 1 ) select a collection of cycles, (Step 
2 ) join some of them together with paths, (Step 3) for each region of the planar dual bounded by 
these cycles and paths, for each cycle G that forms part of the boundary of that region, construct 
the boundary-to-boundary spanner for G-to-G paths in that region. The union of edges from Steps 
1, 2, and 3 form the spanner. 

So far, however, we have not handled weights. Indeed, a perturbation (in which a path P of 
the optimal solution is replaced with a path P' in the spanner) could shift weight from one side of 
the bipartition to the other. We need a way to limit the amount of weight that could shift. This 
is the purpose of Step 1. Note that the original path P and replacement path P' form a cycle G, 
and that weight that could shift is enclosed by G. The goal of Step 1 is to ensure that, for every 
such cycle G derived from such a perturbation, the weight enclosed by G is small compared to the 
cost of G. That way, the total amount of weight shifted can be charged to the cost of the optimal 
solution. 

Step 1 ensures that such cycles’ cost/weight ratios are large by greedily finding a collection of 
mutually noncrossing cycles of small ratio. Once Step 1 has completed, each of the regions bounded 
by the noncrossing cycles contains no cycle with small cost/weight ratio (essentially). 

The fact that cycles found in Step 1 have small cost/weight ratio is used to show that the total 
cost of all those cycles is not much more than the cost of the optimal solution. Each cycle’s cost 
is charged to the weight of some of the faces it encloses. If we ensure that each face is charged at 
most a logarithmic number of times, the total cost of the cycles is at most a log times OPT. 

To ensure logarithmic charging. Step 1 alternates between adding cycles to the spanner and 
removing cycles. When one cycle is enclosed by the other but the two cycles enclose almost the 
same weight, the pair of cycles is designated a splicing pair, and, in an operation called splicing, 
cycles sandwiched between them are removed from the spanner. 

Here is one complication: A region bounded by cycles from Step 1 often is bounded by several 
cycles, i.e. its boundary is disconnected. The boundary-to-boundary spanner (Lemma 2.2) works 
only for a connected boundary. Step 2 therefore uses a technique called PC clustering, due to 
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Bateni, Hajiaghayi, and Marx [6], to add paths joining some of the boundary components. If two 
boundary components remain unconnected after PC clustering, we can assume that some near- 
optimal solution does not connect them. 

Here is another complication: Consider a cycle C associated with a perturbation that replaces 
a path P of the optimal solution with a path P' in the spanner. If C happens to be sandwiched 
between two cycles comprising a splicing pairs, then C might have small cost/weight ratio. This 
happens if C encloses the inner cycle C' of the splicing pair. To make sure no such cycle C is used in 
a perturbation. Step 3, in forming the boundary-to-boundary spanners for the region R containing 
C', distinguishes between paths going clockwise around C' and paths going counterclockwise. This 
is accomplished using a construction from topology, the cyclic double cover. 


3 Preliminaries 

Achieving polynomially bounded weights Garg, Saran and Vazirani [13] observed that if one 
is willing to accept a (6±6:)-bipartition instead of a 6-bipartition, then one can assume polynomially 
bounded weights. This is done by defining new weights w{v) ^ where W is the sum of 

original weights and n is the number of vertices of the input graph. After the transformation, the 
sum of weights is bounded by e~^n. However, due to the truncation, the weight of any vertex may 
be off by ^ with respect to its original weight. Therefore, the weight of any set in a bipartition 
may be off by at most sW. This is allowed by our theorems. We therefore assume henceforth that 
the sum of weights is polynomially bounded. 

Basic definitions We assume that the reader is familiar with the basic concepts of planar graphs 
such as planar embeddings and planar duality. We use G* to denote the planar embedded input 
graph, and we use G to denote its planar dual. The costs of edges in G* are assigned to the 
corresponding edges of G. The weight function can be viewed as an assignment of weights to faces. 
For the remainder of this paper we deal with the dual graph G. Henceforth, unless otherwise stated, 
vertices, faces and edges refer to those of G. 

For each edge e in the edge-set E, we define two oppositely directed darts, one in each orien¬ 
tation. We define rev{’) to be the function that maps each dart to the corresponding dart in the 
opposite direction. 

A non-empty sequence di .. .dk of darts is a walk if the head of di is the tail of for every 
1 < i < k. A walk is said to be a closed walk if the tail of di is the head of dk. 

Let X be a walk in a planar embedded graph, and let P = a X b and Q = c X dhe walks that 
are identical except for their first and last darts. Let a' be the successor of a in P and let F be 
the predecessor of b in P. We say Q forms a crossing configuration with P if the clockwise cyclic 
order of darts whose head is head{a) induces the cycle (a rev {a') c) and the clockwise cyclic order 
of darts whose tail is tail{b) induces the cycle (6 rev{F) d). 

We say a walk P crosses a walk Q if a subwalk of P and a subwalk of Q form a crossing 
configuration. See Figure 1(a). We define a cycle (of darts) to be a non-self-crossing closed walk 
that uses each dart at most once. We omit the modifier “of darts.” Thus for our purposes a cycle 
can use an edge at most twice-once in each direction—and cannot cross itself. 

We can assume (by adding a self-loop if necessary) that G has a face with zero weight whose 
boundary has zero cost. We use foo to denote this face, and we refer to it as the infinite face. 

Let G be the dart multiset of a set of closed walks, and let G* be the multiset of corresponding 
darts in G*. Let / be a face of G. Let P be any foo-to-f path in G* (the input graph). We say / 
is enclosed by G if the parity of |P H G*| — \P H rev{C)\ is odd. See Figure 1(b). 
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(a) Two crossing paths and a cycle of darts (left), a set of mutually noncrossing 
cycles (middle), and a set of cycles that are not mutually noncrossing (right). In 
the middle, the region bounded by three of the cycles is shaded. 



(b) A set C of two cycles (dou¬ 
ble lines). Face a is enclosed 
by C, but face b is not. 

Figure 1. Crossings and enclosure 


Consider a bipartition U UV in G* where the vertex of G* corresponding to /oo belongs to V. 
The duality of cuts and cycles implies that the edges crossing this cut form a set C of cycles in G. 
The weight of the bipartition in G* is the total weight of faces enclosed by C. 

For a cycle G of G, we define wg{G) to be the sum of weights of faces enclosed by G in G. For 
a set S of faces of G, we define ws{G) to be the sum of weights of faces in S enclosed by G, and 
we define rs{G) to be the ratio of the cost of G to ws{G). 

Regions defined by a mutually non-crossing set of cycles Let C be a set of cycles that are 
mutually non-crossing. See Figure 1(a). Assume for the sake of convenience that the boundary of 
/oo is in C. The set C can be represented by a rooted ordered tree T. Every node v corresponds 
to a cycle Gy G C. Ancestry in T is determined by enclosure: Node v is an ancestor of node u in 
T if Gy encloses Gy in G. Thus, /oo is the cycle corresponding to the root of T. 

Every node p of T is associated with a region Ry. Ry is the subgraph of G consisting of vertices, 
edges, and faces enclosed by G^, and not strictly enclosed by Gy for any child u oi v. 

The cycle Gy is called the outer boundary of Ry. For a child u of p, the cycle Gy is called a 
hole of Ry. The weight of a hole is the total weight of faces enclosed by the hole. Together, the 
outer boundary and the holes form the boundary of Ry. We say Ry strietly eontains an edge if in 
addition the edge does not belong to the boundary of Ry. 

We say that a region R eontains a cycle G if i? contains every edge of G, and that R strietly 
eontains G if in addition R strictly contains at least one edge of G. 

Finding low-ratio cycles Let T be a shortest path tree with root r. Let G be a cycle that 
encloses r. We say that G is discovered by T from the inside if, for every p G G, the r-to-v path in 
T is enclosed by G. Rao [26] described a polynomial time technique that finds, for every possible 
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weight w, the minimum-cost cycle enclosing exactly w weight among cycles that go through r and 
are discovered by T from the inside (if such a cycle exists). 

The following lemma implies that Rao’s technique can be used to find a maximally face-enclosing 
cycle with ratio at most some ratio a. 

Lemma 3.1. Let G be a planar embedded graph with edge costs and face weights. Let T be a 
shortest path tree in G, rooted at a vertex r. Let G be a maximally face-enclosing cycle G with 
ratio at most a. If G encloses r then T discovers G from the inside. 

Proof: Assume that G is not discovered by T from the inside. Then there is a vertex v ^ G such 
that the r-to-v path P in T is not enclosed by G. Let P' be a maximal subpath of P consisting only 
of edges that are not enclosed by G. Let x,y be the endpoints of Pf Note that, since G encloses 
r, both X and y are vertices of G. Let Q' be a subpath of G with endpoints x and y such that the 
cycle G' — Q' o P' encloses G. Since P' is a shortest x-to-y path, the cost of G' is at most the cost 
of G. Cycle G' also encloses every face enclosed by G, so the weight of G' is at least the weight of 
G. This contradicts the fact that C is a maximally face-enclosing cycle G with ratio at most a. □ 

To find a maximally face enclosing cycle with ratio at most a, consider the set M of cycles whose 
weight is maximum among all cycles whose ratio is at most a. Let C be a cycle in M enclosing the 
greatest number of faces (there may be faces with zero weight). Note that G is a maximally face 
enclosing cycle with ratio at most a. Let w denote the weight of G. To find G, slightly perturb 
the weight of every face to make it non-zero without significantly changing the total weight of 
any cycle. This can be done by scaling the weights by the number of faces, and adding 1 to the 
weight of every face. Note that this transformation keeps the total weight polynomially bounded. 
For every possible choice of root r of the shortest path tree T, use Rao’s technique to compute, 
for each possible (perturbed) weight x, the minimum cost cycle discovered from the inside, and 
enclosing exactly x perturbed weight. Return the cycle enclosing maximum perturbed weight whose 
unperturbed ratio is at most a. 

Let be the perturbed weight of G. By Lemma 3.1, for a correct choice of r, G is the minimum- 
cost cycle computed for weight w'. By definition of the perturbation, the perturbed weight w' of 
G is greater than that of any other cycle with ratio at most a that encloses fewer faces than G. 
Also, no cycle whose (perturbed) weight is greater than w' has ratio at most a with respect to 
unperturbed weights. Therefore, the procedure described returns the cycle G. 

Edge contractions One step of the algorithm described in this paper uses edge contractions. 
The planar embedding is not relevant in this step. We therefore use a definition of edge contraction 
that does not depend on the embedding. Contracting an edge e results in merging the endpoints 
of e into a single vertex. Any self loops that are created in the process are deleted. 

4 Skeleton construction 

We now begin to describe our procedure for constructing a spanner to approximate a minimum 
cost 6-bipartition. For brevity, we let OPT be the cost of an optimal 6-bipartition. Let W be the 
total weight of all faces in G. Recall that we work directly with the dual graph G of our input 
graph G* and that a solution or 6-bipartition is a set of cycles in G. We assume that the spanner 
algorithm is given a rational number A such that OPT/VF < A < 20PT/VF, since an outer loop 
surrounding the approximation algorithm can try different values of A and return the best solution 
obtained. 
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Figure 2. On the left are non-crossing cycles. On the right is the corresponding tree. Siblings are ordered left to right in 
increasing order of weight enclosed. The cycles are labeled in preorder. 

The algorithm constructs a family C of mutually noncrossing cycles, which includes the outer 
face /oo- We refer to C as the skeleton. As discussed in Section 3, the cycles of C define regions^ 
and define a rooted tree T based on enclosure. For each node v of this tree, we order the children 
of V left to right according to increasing weight of the faces they enclose. See Figure 2. A preorder 
traversal of T that visits siblings in this order defines a total left-to-right ordering on the nodes of 

r. 

The skeleton building algorithm is given as Algorithm 1. As the algorithm progresses, cycles 
and regions are defined with respect to the eurrent set of cycles C. 

Algorithm 1 Skeleton(G, A) 

1: C ^ {boundary of /oo} 

2: ptr ^ null 
3: repeat 

4: if some region R with respect to C strictly contains a cycle C s.t. tr^C) < e~^X then 

5: let C he a maximally face-enclosing such cycle 

6: set F{C) to be the set of faces of R enclosed by C 

7: add C to C 

8 : else 

9: let the cycles of T in preorder he Ci ... Ck 

10: if ptr = null then ptr ^ Ck 

11: else ptr ^ cycle preceding ptr in preorder of T 

12: let Cq be the value of ptr 

13: let p be min such that Cp is an ancestor of Cq and wciCp) < 2wG{Cq) 

14: remove cycles Cp+i,..., Cq-i from C 

15: until ptr = boundary of /oo 
16: return C 


Consider an execution of Line 14 where a sequence of cycles Cp+i,..., Cq-i is removed from C. 
We say that (C^, Cq) is a splieing pair. Cycle Cp is the rootward member of the pair, and Cq is the 
leaf ward member. 

The following lemmas establish several useful properties of the skeleton, including the polyno¬ 
mial running time of the skeleton-building algorithm. 

Lemma 4.1. Let B be the value of ptr at any time t in which Line 12 is executed. The suffix 
of the preorder starting at B does not change after time t. In particular, cycle B will never be 
removed from the skeleton. 
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Proof: The proof is by induction on the number of times line 12 is executed after time t. If this 
number is zero, i.e. if time t is the last execution of line 12, then T can only change when a splice 
with splicing pair (A, B) occurs, for some cycle A. By definition of the algorithm, no new cycles 
can be added to T at time t. Splicing with splicing pair {A,B) does not change the region of any 
cycle to the left of A (exclusive) or to the right of B (inclusive). Therefore, no new cycles can be 
added to the left of A or to the right of B after the splice. In other words, any new cycle C added 
after the splice is added as a descendent of A, but not of B. Since wg{B) > ^wg{A), B always lies 
on the rightmost path in the subtree of T rooted at A. Hence any such new cycle C is added left 
of B in the preorder. 

For the inductive step, by the argument above, the suffix starting at B does not change until 
the next time line 12 is executed. Since ptr moves strictly to the left, the inductive hypothesis 
implies that the suffix remains unchanged for the remainder of the execution. □ 

Lemma 4.2. The skeleton algorithm runs in polynomial time. 

Proof: By Lemma 4.1, each time line 12 is executed, the length of the fixed suffix of the preorder 
of T increases by one. The number of cycles in the skeleton is bounded by the number of faces in 
G, so there are 0{n) splicing steps. Hence, the total number of cycles added to C throughout the 
execution of the skeleton construction algorithm is 0{n^). Since each cycle is found in polynomial 
time, the total running time is polynomial. □ 

Lemma 4.3. The cost of the skeleton is 0(£:“^OPTlog VF). 

Proof: Let / be a face of G, and let G{f) = {G ^ C : f ^ F{G)}. By Lemma 4.1, every cycle A 
in the final skeleton is pointed to by ptr at some time t. Cycle A then participates in a (possibly 
trivial) splice, removing all but one ancestor cycle with weight less than twice wg{A). Afterward, 
for each cycle G added to C, we have F{G) H F{A) = 0. Let Gi,... ,Gk be all cycles of G{f) in 
rootward ordering. The weight of the cycles in the sequence doubles at least once every other cycle. 
We have \G{f)\ = 0(loglF). 

The ratio of each cycle G in the skeleton is bounded by e~^X, so c{G) < e~^XwG{F{G)). Thus, 
the total cost of cycles in the skeleton is at most 

c fee 

< 0{e-^XW\ogW) 

= 0(£“^OPTlogW^). 


□ 

The heavy nesting of low-ratio cycles in the skeleton guarantees the following lemma, which 
will be crucial in proving our spanner construction is correct. 

Lemma 4.4. Suppose that, when the algorithm terminates, a region R contains cycle G ^ C, and 
'^g{C) < £~^X. Then G encloses the heaviest hole of R, and rs{G) > e~^X where S = {faces 
enclosed by outer boundary of R but not by heaviest hole} 

Proof: We say a cycle D weakly crosses GAD crosses G or D — G. The following claim is 
immediate: 

Claim 1: At termination, C contains no cycle that weakly crosses G. 
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Claim 2: Some splice removes a cycle weakly crossing C. 

Proof of Claim 2: At some point the algorithm adds to C a cycle that weakly crosses C or is 
strictly enclosed by C. (If this never happens then C remains available to choose in line 5, and the 
algorithm is not ready to terminate.) Furthermore, by maximality in line 5, before the algorithm 
adds to C any cycle strictly enclosed by C, it must add a cycle that weakly crosses C. By Claim 1, 
such a cycle must have been removed by a splice. □ 

Claim 3: For the last splice removing a cycle D that weakly crosses C, the splicing pair (A, B) 
must satisfy the following condition: 

* each edge of C is enclosed by A and not strictly enclosed by B. 

Proof of Claim 3: Since the splice removes every edge of D is enclosed by A and not 
strictly enclosed by B. After that splice, A ^ C. If C weakly crosses A then by Claim I A must 
be removed later, contradicting the choice of (A,S). Similarly, C does not weakly cross B. This 
implies every edge of C is enclosed by A and is not strictly enclosed by S, proving the claim. 

Claims 2 and 3 imply that there is some splice whose splicing pair satisfies Condition *. Let 
(A^, B^) be the splicing pair of the last such splice. Let to be the time of that splice. 

Let Sq — {faces enclosed by A^ and not by B^}. We claim rsQ{C) > e~^X. If not then after 
time to the cycle C would still be available to add to C, so some cycle weakly crossing C would 
have been added, contradicting either Claim I or the definition of to. Because rs{C) > ra{C) and 
B^ strictly contains no edge of C, it follows that C encloses B^. By Lemma 4.1, cycle B^ is never 
subsequently removed from C. Also, by Claims I and 3 no cycle weakly crossing C is subsequently 
added to C. Therefore, for every t > to, at time t there is a unique pair A^, B^ of cycles in C such 
that A^ encloses (7, and C encloses and A^ is the parent of B^ in the tree T of cycles of C. 

Since B^ remains in C, we infer B^ encloses B^ for all t > to- We show by induction that at 
each time t > to, cycle A^ encloses Ab There are two mutually exclusive cases. First, suppose that 
just after time t a cycle A^+^ enclosing C is added to C and becomes the parent of B^. Then, since 
A^ is still in C, A^ encloses A^+^ and A^ encloses A^+^. 

Second, we show that no splice at time t + I can remove Ab Assume for a contradiction that 
some splice did remove A^, and let (A, S) be the splicing pair. If B were not a descendant of A^ 
in T then the fact that the splice removes A^ would mean it would remove all descendants of A^, 
including a contradiction. If B crosses C then by Claim 3 the splice or a later one satisfies 
Property *, contradicting the choice of to- Therefore S is a proper descendant of A^ that does not 
cross C. B cannot enclose C else it would be a parent of contradicting the fact that A^, B^ are 
a parent-child pair. Therefore B strictly encloses no edge of C. But then (A, B) satisfy Property *, 
a contradiction. 

Let T be the time the algorithm terminates. We have shown that A^ encloses A^ and B^ 
encloses B^. Let S = {faces enclosed by A ^ and not by B^}. It follows that wg{B^) > wg{B^) 
and wg{^) < wg{A^) and S So> Because wg{A^) < 2wg{B^)., we have wg{^) < 2wg{B^)^ 
so B^ is the heaviest hole of the region whose outer boundary is A^. Because S C So and 
rso{C) > e“^A, it follows that rs{C) > e~^\. This completes the proof. 


5 Shortcuts 

Consider an optimal solution O. Let K he a cycle of O that crosses the skeleton. A path decomposi¬ 
tion of K is a decomposition of K into paths • • • such that the number of paths is minimized 

and no path pi crosses the skeleton. Because the number of paths is minimized, the endpoints of 
the paths occur only on skeleton cycles that are crossed by K. 
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Figure 3. (a) Path r of a path decomposition enters and leaves a region at vertices u and v. Replacing r with shortcut q 
will only affect the enclosure of faces outside the largest hole of the region. Replacing r with shortcut p will affect faces 
within the hole as well, (b) The cyclic double cover of the region. Lifts of r and q have the same endpoints in the cyclic 
double cover, but no lift of p has the same two endpoints. Therefore, p would not be a considered a shortcut for r within 
the cyclic double cover. 


Let be a member of a path decomposition as defined above. Let u and v be the endpoints 
of p and let R be the region of the skeleton that contains p. In order to build the spanner, we 
wish to replace p by a (1 + 6:)-approximately shortest 'u-path p' in R. For concision, we call the 
approximately shortest path p' a shortcut The spanner contains a set of edges that carry shortcuts 
for all paths of all path decompositions. These edges are chosen so their total cost is sufficiently 
small and replacing each path by its shortcut does not perturb the weight enclosed by the solution 
by more than a small amount. 

If R contains no holes, then it is relatively easy to compute shortcuts within R. Let C be the 
outer boundary of R. The spanner algorithm computes a subgraph A of G within R using the 
boundary-to-boundary spanner of Lemma 2.2. 

As we show later in Lemma 8.3, replacing paths through regions without holes by their shortcuts 
does not change the weight of the solution by very much. Essentially, the faces that change sides 
when a path p is replaced by its shortcut p' are exactly those enclosed by p o rev{p'). The cycle 
enclosing these same faces must have low weight or there would exist a low ratio cycle within R 
violating Lemma 4.4. 

Unfortunately, this argument does not hold if R contains holes. First, there may be paths in 
a path decomposition that start and end on different boundary components of R. A boundary-to- 
boundary spanner does not contain shortcuts for such paths. And even if shortcuts are available, po 
rev{p') could enclose a high weight hole, and all of the faces within that hole could change sides 
after replacing p with p'. See Figure 5 (a). Lemma 4.4 does not imply anything about a cycle’s 
weight if it encloses the largest hole of a region, so there would be no limit to how much weight 
could change sides. In the next section, we describe how to address both of these issues for regions 
with holes. We begin by describing a topological construct called the cyclic double cover used to 
address the later issue. A procedure called PC-clustering will be used to address the former. 


6 The cyclic double cover 

In this section, we describe a tool called the cyclic double cover that is used by our spanner 
algorithm. Thanks to this tool, the spanner will carry shortcuts that do not force the largest 
hole in each region to change sides. The cyclic double cover as described here was originally used 
by Erickson [9] and by Fox [11] to find short topologically interesting cycles in surface embedded 
graphs. Our presentation of the cyclic double cover is based closely on theirs. 







12 


A Polynomial-time Bicriteria Approximation Scheme for Planar Bisection 


Let i? be a region with outer boundary Cq and largest hole Ci. Let L be an arbitrary path 
from Co to Ci. Cut along L to create a new region R' where Cq, and two copies of L (L+ and 
L~) form the outer boundary. Let (i?^0) and (i?', 1) be two distinct copies of R'. For any vertex 
V in i?, let (^’,0) denote the copy of v in and let 1) denote the copy in (i?', 1). Finally, 

let (L^,0) denote the copies of in and (L^,l) denote the copies in The cyclic 

double cover R^ is the planar graph resulting from identifying (L+, 0) and (L“, 1) to a single path 
(L,0) and identifying (L+, 1) and (L“,0) to (L, 1). Every hole of R appears twice as a hole in R^ 
except for Ci. The edges of Cq appear twice along the outer face of R^ and the edges of Ci appear 
twice along a single hole in R^. See Figure 5 (b). 

The cyclic double cover R^ has an equivalent combinatorial definition. Each vertex v from R 
has two copies ('r’,0) and 1) in R^. Edge uv has two copies (?i, 0) and 1) if uv does 

not enter L from the left. Otherwise, the copies of uv are (i/,0)('?;, 1) and (?i, !)('?;, 0). Edge copies 
in R^ retain their costs from R. The projection of any vertex, edge, or walk in R? is the natural 
map to R that occurs by dropping the 0 or 1 from the vertex and edge tuples. We say a vertex, 
edge, or walk p in R lifts to p' if p is the projection of p'. The outer boundary of R^ projects to 
two copies of Cq and one boundary of R^ projects to two copies of Ci. Otherwise, every face or 
boundary of R^ projects to a face or hole in R. Call a boundary in R^ a hole if it projects to one 
or two copies of a hole in R. For a walk p in i?, let x^(p) = 0 if p crosses L an even number of 
times. Otherwise, let x‘^{p) = 1. We immediately get the following lemmas. 

Lemma 6.1. Let p be a walk in R from vertex u to vertex v. Walk p is the projection of a unique 
walk p^ in R? from vertex (u,0) to vertex (v,x^(p)). 

Lemma 6.2. Let C' be a closed walk in R^. C' projects to a unique closed walk C in R that does 
not enclose Ci. 

Our algorithm builds a spanner by computing the cyclic double cover for each region R of the 
skeleton that contains a hole. The construction of the double covers can be done in linear time. It 
then computes a set of edges within each double cover that carry shortcuts for the lifts of paths 
that may appear in path decomposition. See Figure 5 (b). The projections of these edges back 
into the original regions will be added to the spanner. In Lemma 8.4, we will show that we can 
replace the paths of path decompositions by the projections of these shortcuts without changing 
the weight of the solution by more than a small amount. 

7 PC-clustering 

Our spanner algorithm needs to find edges that carry shortcuts within each cyclic double cover. 
However, the 

boundary-to-boundary spanner of Lemma 2.2 will not find edges to carry shortcuts between distinct 
boundaries of a region. In order to use the algorithm, we augment the skeleton’s edges within each 
double cover using the PC-clustering algorithm of Bateni, Hajiaghayi, and Marx [6]. 

Let K be a cycle in the optimal solution that crosses the skeleton, and let p be a path of iF’s 
path decomposition where p lies in a region R with a hole. The PC-clustering algorithm adds a 
relatively cheap set of edges to the boundary of the cyclic double cover R^. These edges are chosen 
so that, in general^ both of p’s endpoints lie on the same boundary component after running PC- 
clustering. If p’s endpoints are still on different components, then at least one of the components 
must be very cheap. The edges of that component can be added to the optimal solution without 
substantially increasing its cost, and K can be modified to avoid crossing the skeleton cycles in 
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that component. Path p is no longer in a path decomposition, and it is no longer necessary for the 
spanner to hold a shortcut between p’s endpoints. Formally, the PC-clustering algorithm can be 
described as follows. 

Lemma 7.1 (Bateni et al. [6]). Let G{V,E) be a graph with non-negative edge costs c(e) and 
face potentials There exists a polynomial time algorithm to hnd a subgraph Z such that 

1. the total cost of Z is at most 2 J2vev 

2. for any subgraph H of G, there is a set U of vertices such that 

(a) is at most the cost of H and 

(b) if two vertices vi,V 2 ^ U are connected by H, then they are in the same component of 

Z. 

For each region R containing a hole, our algorithm does the following: It contracts the lifts of 
every skeleton edge in to get the graph R^. For any vertex v in K^, let c{v) be the total cost of 
edges contracted to create v (implicitly, if v appears in R^ as well, then c{v) = 0). For each v in IP, 
the algorithm assigns a potential (l){v) — The algorithm them applies the PC-clustering 

procedure of Lemma 7.1 to R^ to get the set of edges Z. 

Let the well-connected cover graph be the set of boundary in R? unioned with Z. The edges 
in the well-connected cover graph are the well-connected cover edges. The projections of the well- 
connected cover edges will be used in the spanner. 

Lemma 7.2. The total cost of all well-connected cover edges is 0{e~^OPTlogW). 

Proof: By Lemma 4.3, the total cost of all cycles in the skeleton is 0(6:“^0PT log W). Every cycle 
in the skeleton appears on the boundary of at most two regions. For each boundary of a region, each 
edge of the boundary appears at twice in that region’s cyclic double cover. Therefore, the sum of 
vertex potentials used for PC-clustering across all cyclic double covers is at most 0{s~^0PT log W). 
The lemma follows from the first property of PC-clustering’s output as defined in Lemma 7.1. □ 

We argue there exists a near-optimal solution such that shortcuts do not start and end on 
different components of the well-connected cover graphs. 

Lemma 7.3. There exists a solution with cost at most (1 -h 2£:)OPT enclosing exactly bW weight 
such that for each cycle K in the solution, 

1. either K does not cross the skeleton 

2. or K has a path decomposition po,pi,... such that a lift of each path pi in a region with a 
hole has both endpoints on the same component of that region’s well-connected cover graph. 

Proof: Consider the optimal solution O. For each region R with a hole, let Or be the subset of 
edges from O that lie strictly within the region R. Consider a lift of Or to R^, and let Or be the 
edges that remain after performing contractions to get R^. Let Ur be the set of vertices in IP that 
are guaranteed to exist for Or by the second property in Lemma 7.1. 

Each vertex of Ur is the result of contracting zero or more lifted skeleton cycles in R^. Let U be 
the set of all skeleton cycles where for each cycle O ^U, the contraction of a lift of O lies in some 
Ur. The cycles in U are mutually non-crossing, so they partition the faces of G into U-regions. 
For each ZY-region R', let Or> be the boundary of faces in R' enclosed by O (therefore, the holes 
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of i?' are not enclosed by Let O' be the union of cycles over all Oi^/. Solution O' encloses 

the same set of faces of G as O. No cycle of O' crosses a member of U. Edges strictly internal 
to some Z//-region are used exactly once and only if they are used in O. Finally, each cycle C 
may contribute up to two copies of some of its edges to O', because C lies on the boundary of two 
Z//-regions. 

The cost difference between O' and O is at most twice the cost of cycles in U. For a single 
region i?, the total potential of vertices in JJr is at most the cost of Or. By definition of vertex 
potentials, the cycles of U that contribute to set JJr have total cost at most £c{Or). All sets Or 
are edge-disjoint, so twice the total cost of all cycles in U is at most 2£:OPT. 

Now, consider any cycle K in O' that crosses the skeleton, and let • • • be a path decom¬ 

position for K. The endpoints of each path pi lie on skeleton cycles crossed by K. For any pi lying 
in a region R with a hole, let p[ be a lift of path pi to R^. Let be the path that results from 
p' after contracting edges to make R^. Finally, let vi and V 2 be the endpoints of Neither vi 
nor V 2 lie on a member of Ur as defined above, because the cycles of O' do not cross any members 
of U. Further, every edge of p' that is not contracted is a member of Or. Therefore, Lemma 7.1 
guarantees vi and V 2 lie in the same component of PC-clustering’s output. Each vertex of W is a 
connected component of the boundary of R? so the endpoints of p' lie on the same component of 
the well-connected cover graph as well. □ 

7.1 Finding shortcuts 

For each cyclic double cover i?^, for each well-connected cover component, our spanner algorithm 
computes edges carrying shortcuts between every pair of vertices on the component. It does so using 
the following extension of the boundary-to-boundary spanner of Lemma 2.2. The projection of these 
edges is added to our spanner. 

Lemma 7.4. Let G be a planar graph with non-negative edge costs c(-). Let A be a component of 
G. Then for any e > 0, there is an 0{n log n) time algorithm to compute a subgraph H of G where 
c{H) = 0(6:“^c(A)) and for any pair of vertices u and v on A, dist/f (?i, 'L’) < (1 + 6:)distG(?i, v). 

Proof: The edges and vertices of A partition the plane into one or more components. The algorithm 
cuts the planar graph G along A, separating these components so that the edges of A lie on their 
boundary. For each boundary cycle G of the cut open graph, the algorithm runs the boundary- 
to-boundary spanner procedure of Lemma 2.2 to create a subgraph He of G of cost at most 
0{£~^c{G)) in 0{nclognc) time where nc is the number of vertices in (7’s component of the 
cut open planar graph. Subgraph H is the union of all such subgraphs He- Each edge of A 
appears on boundary twice, so the total cost of all subgraphs He is at most 0(6:“^c(A)). The 
cut open surface has 0{n) vertices total, so the running time for the procedure is O(nlogn) total. 
Finally, the boundary-to-boundary shortcut algorithm guarantees that for any pair of vertices u 
and 'L’ on A where the shortest path does not cross any cycle G of the cut open surface, we have 
distniu^v) < (1 + 6:)disto'(i/,'L’). This proves the lemma since a shortest path that does cross a 
cycle G of the cut open surface is the concatenation of shortest paths that do not cross any such 
cycle. □ 

Lemma 7.5. The total cost of all edges carrying shortcuts is 0(6:“^0PT log VF). 

Proof: By Lemma 4.3, the total cost of boundaries for regions without a hole is 0(£:“^0PT log W). 
By Lemma 7.2, the total cost of all subgraphs A used in Lemma 7.4 for regions with holes is 
0(6:“^OPTlog VF). The current lemma follows from Lemmas 2.2 and 7.4. □ 
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8 The spanner 

In this section, we describe the final spanner construction for our algorithm and prove the con¬ 
struction follows the spanner properties. The spanner construction is summarized as Algorithm 2. 
The cheapest cycle other than the outer boundary enclosing a particular hole of a region R can be 


Algorithm 2 Spanner(G) 

1: 5^0 ; C ^ Skeleton(G, A) 

2: add to S the edges of C 

3: for each region R without a hole, add to S the edges carrying shortcuts in R 

4: for each region R with a hole, add to S the projection of the well-connected cover edges from R^ 

5: for each region R with a hole, add to S the projection of the edges carrying shortcuts in R^ 

6: for each region R with a hole, add to S the edges on the cheapest cycle in R that encloses the 
heaviest hole of R other than the outer boundary of R 


computed in polynomial time using several instantiations of any polynomial time minimum sR-cut 
algorithm [15]. 

To prove that our algorithm computes a spanner, we will iteratively replace cycles and their 
paths in a near-optimal solution with ones that lie in the spanner. The following lemmas will help 
us bound the total change in weight and cost from performing these operations. 

Lemma 8.1. Let O be a set of cycles. Let K ^ O be contained by region R, and let rciK) > e~^X. 
Removing K from O changes the enclosed weight of O by at most ec{K)/\ and does not increase 
the cost of O. 

Lemma 8.2. Let O be a set of cycles. Let K ^ O be strictly contained by region R, and let 
'^g{K) < e~^\. Cycle K encloses the largest hole of R. Further, replacing K with the shortest 
cycle other than the outer boundary of R that encloses the largest hole of R changes the weight 
enclosed by O by at most ec{K)/\ and does not increase the cost of O. 

Proof: By Lemma 4.4, we know K must enclose the largest hole of R. Let A be the outer boundary 
of R and B be the largest hole of R. Let K' be the shortest cycle other than A enclosing B in R. 
Replacing A by a possibly cheaper cycle K' cannot increase the cost of O. Let Aw{K) be the total 
weight of faces that move from being enclosed to not enclosed and vice versa after the replacement. 
Let Aw{K) (respectively Aw{K')) be the total weight of faces that are enclosed by K [K') but 
not enclosed by B. Suppose Aw{K) > Aw{K'). By Lemma 4.4, we have c{K)/Aw{K) > e~^\. 
In particular, Aw{K) < Aw{K) < ec{K)/X. If Aw{K') > Aw{K), then either Aw{K') = 0 or K' 
is not in C. Either way, we have c{K')/Aw{K') > e~^X, and Aw{K) < Aw{K') < ec{K')/X < 
ec{K)/X. □ 

Lemma 8.3. Let O be a set of cycles. Let K ^ O cross one or more cycles of the skeleton. Let 
p be a path of K^s path decomposition lying in region R such that R has no holes and p has at 
least one edge disjoint from the spanner. Let u and v be the endpoints of p on the boundary of R. 
Finally, let p' be a shortcut in the spanner between u and v in R. Replacing p by p' changes the 
weight enclosed by O by at most 3ec{p)/X and increases the cost of O by at most ec{p). 

Proof: The bound on the cost change follows from the fact that shortcut p' is a (l-h6:)-approximate 
shortest path. Let e be the edge of p strictly internal to R. Let C — po rev{p'). Let C' be the cycle 
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using darts of C that encloses and connects the boundary of faces of G enclosed by C. Cycle C' is 
strictly enclosed in i?, because it contains e. Also, c{C') < c{C). The faces of G enclosed by G' are 
exactly the faces that switch sides when replacing p hy p'. Therefore, the replacement will change 
the weight enclosed by O by exactly wg{G'). By Lemma 4.4, we have rciG') > e~^X. Recall, a 
shortcut is a (1 + £:)-approximate shortest path. Therefore, 

.-1A < = c(p) + c{p') ^ (2 + g)c(p) 

waiCO wg{C0 - wg(C') • 

In particular, wciC') < e(2 -\- e)c{p)/\ < 3ec{p)/X. □ 

Lemma 8.4. Let O be a set of cycles. Let K ^ O cross one or more cycles of the skeleton. Let 
p be a path of K^s path decomposition lying in region R such that R has at least one hole and p 
has at least one edge disjoint from the spanner. Let p' be a lift of p from R to R^ with endpoints 
u and V. Finally, let p" be a shortcut in the spanner between u and v in R^. Replacing p by the 
projection of p" changes the weight enclosed by O by at most 3ec{p)/X and increases the cost of O 
by at most ec{p). 

Proof: Again, the bound on cost change is immediate. Let G be the projection of p' o rev{p''). By 
Lemma 6.2, (7 does not enclose the largest hole of R. The rest of the proof is identical to that of 
Lemma 8.3. □ 

We finally prove that our algorithm constructs a spanner. 

Lemma 8.5. The output of Algorithm 2 contains a solution with cost at most (1 + 46:)OPT 
enclosing b'W weight with F e [b — be, b-\-be\. 

Proof: Let O be the near-optimal solution given by 

Lemma 7.3. We will modify O to create a solution O' that lies in the spanner. 

Consider any cycle K in O. If it matches the hypothesis of Lemma 8.1, then it is simply removed 
from O. If it matches the hypothesis of Lemma 8.2, then it is replaced by the shortest cycle other 
than the outer boundary of its region enclosing the largest hole of its region. Line 6 of Algorithm 2 
guarantees this cycle exists in the spanner. 

In the final case, we know K has a path decomposition p{),pi, .... We replace each pi in the 
decomposition as follows. If pi has no edges disjoint from the spanner, then it remains as is. If pi 
has an edge strictly internal to a region without holes, then we replace pi with a (I+6:)-approximate 
shortest path between its endpoints. Line 3 guarantees this path exists in the spanner. Finally, if pi 
has an edge strictly internal to a region R with a hole, then p'-, lift to R^, has both endpoints on 
the same component of i?’s well-connected cover graph. The output of the algorithm in Lemma 7.4 
contains an (I -h £)-approximate shortest path p'' in the cyclic double cover between the endpoints 
of p'. Lines 4 and 5 of Algorithm 2 guarantee the projection of p'j exists in the spanner. Path pi 
is replaced by the projection of p'l. Any remaining cycles or decomposition paths of O exist in the 
skeleton itself, which is also in the spanner by Line 2 of Algorithm 2. 

Solution O' is formed by performing the replacements described above. Let p be any cycle or 
path replaced above. By Lemmas 8.2, 8.3, and 8.4, after replacing p, the cost of O increases by at 
most ec{p)/X, and the weight of faces enclosed by O changes by at most 3ec{p)/X. By summing 
over all cycles K, we see the total cost increases by at most e{l -\- 26:)OPT < 2^0PT (for <l/2) 
and the weight changes by at most 3^(1 + 2£:)OPT/A < 6eW. □ 

Lemma 8.6. The output of Algorithm 2 has cost at most 0(6:“^0PT log VF). 
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Proof: The shortest cycles enclosing largest holes within their region have total cost at most that 
of the skeleton itself. The rest of the spanner’s components have their costs bounded in Lemmas 4.3, 
7.2, and 7.5. □ 

Final remarks We have given a bicriteria approximation scheme, but there is much room for 
improvement. Is there an efficient PTAS, one whose running time is a polynomial with degree 
independent of e? Is there a PTAS that does not approximate the balance b? Is there perhaps even 
a polynomial-time algorithm for bisection in planar graphs when the weights are small integers? 

The research described in this paper began at a meeting at Brown University’s Institute for 
Computational and Experimental Research in Mathematics. The authors thank Daniel Marx who 
described to the second author some results on fixed-parameter tractability of graph separation, 
leading to some discussion of the complexity in planar graphs and MohammadTaghi Hajiaghayi 
who asked whether there is an approximation scheme for bisection in planar graphs. Hajiaghayi 
had some initial thoughts due to its similarity to multiway cut for which there is a PTAS on planar 
graphs. 
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